Virial corrections to simulations of heavy ion reactions 
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Within QMD simulations we demonstrate the effect of virial corrections on heavy ion reactions. 
Unlike in standard codes, the binary collisions are treated as non-local so that the contribution of 
the collision flux to the reaction dynamics is covered. A comparison with standard QMD simulations 
shows that the virial corrections lead to a broader proton distribution bringing theoretical spectra 
closer towards experimental values. Complementary BUU simulations reveal that the non-locality 
enhances the collision rate in the early stage of the reaction. It suggests that the broader distribution 
appears due to an enhanced pre-equilibrium emission of particles. 



The Boltzmann equation including the Pauli blocking 
(the BUU equation Q) and the closely related method 
of quantum molecular dynamics [^]|| (QMD) are exten- 
sively used to interpret experimental data from heavy ion 
reactions. Due to their quasi-classical character, they of- 
fer a transparent picture of the internal dynamics of re- 
actions and allow one to link observed particle spectra 
with individual stages of reactions. 

An ambition to cover the heavy ion reactions within 
experimental errors has been recently cooled down by a 
failure of BUU simulations to describe the energy and 
angular distribution of neutrons and protons in low and 
mid energy domain Indeed, the Boltzmann equa- 

tion is not the full story. As noticed in numerical studies 
of hard sphere cascade by Halbert (?]] and more general 
by Malfliet it is disturbing that all dynamical mod- 
els rely more or less on the use of the space- and time- 
local approximation of binary collisions inherited from 
the Boltzmann equation. This approximation neglects a 
contribution of the collision flux to the compressibility 
and the share viscosity which control the hydrodynamic 
motion during the reaction. To include the collision flux 
and other virial corrections, the non-local character of 
binary collisions has to be accounted for. Malfliet also 
demonstrated that non-local collisions can be easily in- 
corporated into simulation BUU codes. 

In absence of a first principle theory, Malfliet in his pio- 
neering study, and more recently Kortemeyer, Daffin and 
Bauer H , had to use classical hard-sphere-like non-local 



collisions which naturally do not result in a full quan- 
titative agreement with experimental data. This ad hoc 
approximation reflects a gap in former quasi-classical the- 
ories of quantum transport: authors either cared about 
non-local collisions leaving aside quasiparticle features or 
vice versa. Moreover, quantum theories of binary colli- 
sions treated non-local collisions via gradient contribu- 
tions to the scattering integral which are numer- 
ically inconvenient and thus have never been employed 
in demanding simulations of heavy ion reactions. Re- 
cent theoretical studies have filled this gap. Danielewicz 
and Pratt jl2| pointed out that the collision delay can 
be used as a convenient tool to describe the virial correc- 
tions to the equation of state for the gas of quasiparticles. 
Although their discussion is limited to the equilibrium, 
it marks a way how to introduce virial corrections also 
to dynamical processes. This task was approached in 
|l3| , where a quasi-classical kinetic equation was derived 
by a systematic quasi-classical approximation of non- 
equilibrium Green's function in the Galitskii-Feynman 
approximation. The derived kinetic equation is suitable 
for simulations having the quasiparticle renormalizations 
in the standard form of Landau's theory and non-local 
collisions reminding classical hard spheres. 

In this letter we follow the line initiated by Malfliet, 
however, with more advanced theoretical and numerical 
tools. We start from the kinetic equation from jl3| and 
evaluate non-local collisions from the Paris potential solv- 
ing the two-particle T-matrix. These non-local correc- 
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tions are then incorporated into the QMD simulation. 
Results show that the non-local corrections bring the en- 
ergy distribution of protons closer to experimental values. 
On a complementary BUU simulation we enlighten the 
microscopic mechanism leading to this improvement. 

The scattering integral of the non-local kinetic equa- 
tion derived in jl3| corresponds to a following picture of 
a collision. Assume that two particles, a and b, of ini- 
tial momenta k and p start to collide at time instant t 
being at coordinates r a and r^. Due to a finite range of 
the interaction, at the beginning of collision particles are 
displaced by r& — r a = A bc . The collision has a finite 
duration A t , i.e., it ends at t + A t . During the colli- 
sion, both particles move so that their end coordinates 
differ from those at the beginning, r' a — r a = A a and 
r' b — Tf, = Af,. The particle a transfers a momentum 
q to the particle b, therefore their relative momentum 
changes from k = ^(k — p) to k' = i(fc — p) — q. Their 
sum momentum is modified by an external field acting 
on the colliding particles during the collision going from 
K = k + ptoK' = k+ p + Ak- The same field changes 
the sum energy of colliding particles from E = e a + (b to 
E' = e' a + e' b = e a + e b + A E . 

The values of A's are given by derivatives of the total 
scattering phase shift <f> = Im lnTfl(f2, k,p, q, t, r) of the 
two-particle T-matrix Tr, 
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Note that energy Q enters as an independent quantity so 
that one needs to know the scattering phase shift out of 
the energy shell. The on-shell energy, — E' . is sub- 
stituted after derivatives are taken. Figure [l] illustrates 
the nonlocal concept derived in jl3| . We like to point out 
that this concept leads to a continuous trajectory in the 
kinetic picture replacing real potential scattering. Con- 
sequently the energy, momentum, density and angular 
momentum conservation are conserved including second 
order quantum virial corrections fl3|| . 

It is our intention to incorporate these features of col- 
lisions into the QMD and BUU simulation codes. The 
selfconsistent evaluation of all A's for all collisions would 
be too demanding. We employ two kinds of additional 
approximations. First, following approximations used 
within the BUU equation, we neglect the medium effect 
on binary collision, i.e., use the well known free-space T- 
matrix. Second, we rearrange the scattering integral into 
an instant but non-local form. This instant form paral- 
lels hard-sphcre-like collisions what allow us to employ 



computational methods developed within the theory of 
gases O] similarly as it has been done in Q . 

In the instant approximation we let particles to make a 
sudden jump at time t from r a and rb to effective final co- 
ordinates f a and rv These effective coordinates and mo- 
menta R and K are selected so that at time t+A t particles 
arrive at the correct coordinates, r' a and r' b , with the cor- 
rect momenta, k' and K' . Accordingly, in the asymptotic 
region, after t + A t , there is no distinction between the 
non-instant and instant pictures, which is shown as solid 
line in figure |^. This asymptotic condition is naturally 
met if one extrapolates the out-going trajectories from 
known coordinates and momenta at t + A t back to the 
time t. Doing so one finds that the effective coordinates 
read 

, k — q . . kj q , . , 

r„ = r'„ ^A 4 = r a + A a ^A t , (2) 
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The change of the relative momentum is insensitive to 
the instant approximation, k = k' = k — q. The sum 
momentum and energy, however, get simplified because 
during the instant process the mean-field has no time 
to pass any momentum and energy to the colliding pair. 
Accordingly, K = K and E = E. In other words, the 
momentum and energy gains are naturally covered by 
the effect of the mean field on particles during the time 
interval (i, t+A t ) which, in the instant picture, is already 
covered by a free motion. Similarly, in agreement with 
the continuity of the center of mass motion, one finds 
that f a + h = r a + n. 

When incorporating the displacements into the QMD 
simulation code, we have to face the fact that two par- 
ticles are selected for a collision if they meet at the 
point of closest approach. This distance is different from 
the distance A be required from the equivalent scattering 
scenario presented in figure |^ as solid line. We consider 
now the time required to travel from A^e to the dis- 
tance of closest approach A t = ^jkA^c in analogy to 
|l5| . Within this scenario we are allowed to jump at the 
point of closest approach to the final asymptotics ®) and 
(||)whith the additional distance the particle travel dur- 
ing A t . The effective final coordinates thus have to be 
evaluated as 



R a + Rb . f 

Ta.b = ^ T A , 



with the effective displacement 
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Since the center of mass does not jump in the collision, 
the final displacement can be also written in an alterna- 
tive way, A f = iA bc + A b - ^(A t - A t ). The non- 
local corrections are thus performed as follows. When 
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the collision is selected, we evaluate A f from (||) and ([!]) , 
redisplay particles into r a and fb and continue with the 
simulation. 

At this point it is possible to establish a connection 
of the present theory to the hard-sphere-like corrections 
used by Malfliet ||| and Kortemeyer, Daffin and Bauer 
H . For hard spheres of the diameter d, the phase shift 
has a classical limit (j> = 7r — \q\d which gives A a ^ = 
and A f = A bc = -^d. The displacement thus has the 
same amplitude d for all binary collisions and points in 
the direction of the transferred momentum, as it is known 
from the Enskog equation Jig] . In the present theory, an 
amplitude of displacement^) depends on the selected 
channel, A l (k,p, q), and the direction does not coincide 
with the transferred momentum, 

The second term is perpendicular to the transferred mo- 
mentum and stays in the collision plane. The notation 
of coefficients iJi.2 is identical to the corresponding in- 
terface subroutines which can be obtained from authors. 

In order to investigate the effect of non-local shifts on 
realistic simulations of a heavy ion reaction, we have eval- 
uated A f from the two-particle scattering T-matrix T R 
in the Bethe-Goldstone approximation jl7]|ll| using the 
separable Paris potential The comparison of the 

shifts calculated for different potentials concerning par- 
tial wave coupling up to D-waves can be found in |L9| |. 
We have incorporated these shifts into a QMD code for 
the central collision of 129 Xe^ 119 Sn at 50 MeV/A. Fig- 
ure [|a shows the exclusive proton spectra subtracting 
the protons bound in clusters. This procedure is per- 
formed within a spanning tree model which is known to 
describe a production of light charged cluster in a reason- 
able agreement with the experimental data, Figs. |3^z>-f . 
Within the local approximation, however, the remaining 
distribution of high-energy protons is too low to meet 
the experimental values. As one can see, the inclusion 
of non-local collisions corrects this shortage of the QMD 
simulation. As demonstrated in Fig. ^, productions of 
light clusters are rather insensitive to the non-local cor- 
rections. This also shows that the improvement of the 
proton production is not on cost of worse results in other 
spectra. 

A microscopic mechanism leading to the increase in 
the high-energy part of the particle spectrum can be 
traced down to an enhancement of the number of col- 
lisions at the pre-equilibrium stage of the heavy ion re- 
action demonstrated in Fig. ^ for the BUU simulation 
of the same reaction. This enhancement gives rise to 
an immediate proton production which itself translates 
into a high energetic spectra. In other words, the strong 
production of the high-energy protons follows from the 
pre-equilibrium emission of particles. The BUU simula- 
tion also shows that non-local corrections are important 



namely in the early stage of reaction well before most of 
light clusters are formed. It explains why the production 
of protons is affected while the formation of light clusters 
is nearly untouched by the non-local corrections. 

In summary, as documented by the improvement of the 
high-energy proton production, the non-local treatment 
of the binary collisions brings a desirable contribution 
to the dynamics of heavy ion reactions. According to an 
experience from the theory of gases, one can also expect a 
vital role of non-localities in the search for the equation 
of state of the nuclear matter. It is encouraging that 
the non-local corrections are easily incorporated into the 
BUU and QMD simulation codes and do not increase 
computational time. 
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FIG. 2. A nonlocal binary collision (solid line) together 

with the szenario of sudden jump at the closest approach. 
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FIG. 3. The particle spectra for central collision of 129 Xe— > 119 Sn at 50 MeV/A with and without non-local corrections. The 
data are extracted from recent INDRA experiments |}(| . The non-local corrections bring the spectrum of the protons towards 
the experimental values leaving the clusters almost unchanged. 
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FIG. 4. The number of collisions per time with and without non-local collisions within a BUU simulation of the same reaction 
as in figure (^). 
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